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A Numerical Study of Evanescent Fields in Backward- Wave Slabs 
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Numerical study of evanescent fields in an isotropic backward-wave (BW) (double negative or 
"left-handed") slab is performed with the FDTD method. This system is expected to be able to 
restore all spatial Fourier components of the spectrum of a planar source, including evanescent 
fields, realizing a "superlens". The excitation of surface modes on the interfaces of the slab, which 
is the key process responsible for the sub-wavelength focusing, is numerically confirmed and the 
time-domain field behavior is studied. In particular, a numerical verification of the amplification of 
evanescent modes by an isotropic BW slab with e r — fi r = — 1 is provided. 

PACS numbers: 02.70.Bf, 41.20.Jb, 77.22.Ch, 77.84.Lf 



The resolution wave limit of any optical device is well 
known: it is impossible to resolve details smaller than 
half wavelength. The physical reason of this limitation 
comes from the fact that evanescent waves in the Fourier 
spectrum of an object exponentially decay in the di- 
rection from the object plane. The decay factor reads 
a = y/kf — k 2 , where kg is the wavenumber in free space, 
and k t is the wavenumber of Fourier components in the 
object plane. The faster the field varies in the object 
plane, the faster it decays in the direction normal to the 
object plane. However, it was recently found that a very 
special kind of lens made of a material with negative 
relative parameters e r = — l,/x r = —1 (at a certain fre- 
quency) can "amplify" the evanescent part of the spec- 
trum, thus opening a way to realize a "superlens" with 
sub- wavelength resolution 0, 0, . 

The reason for this counterintuitive behavior is the fact 
that for a fixed frequency, an interface between free space 
and a backward- wave medium with e r = — X,ji r = —1 
supports surface waves with arbitrary propagation con- 
stants along the interface 0]. Indeed, the eigenvalue 
equation for surface modes (surface polaritons) on an in- 
terface between two isotropic media with parameters 
and ^x,2 reads 
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where indices 1,2 correspond to the two media. If 
€2 = — ei and fX2 = —fii, both ^ and are identically 
satisfied for all propagation constants kt along the in- 
terface. This means that any evanescent incident plane 
wave is exactly in phase with one of the eigenmodes of 
the surface wave spectrum. If the interface is infinite 
in space, the amplitude of the excited surface wave be- 
comes infinite. The amplification of evanescent fields in 
backward-wave slabs utilizes this resonant excitation of 
waveguide modes with large propagation constants. The 
incident field excites an eigenmode of the slab that is 



formed by two exponentially decaying field components 
inside the slab. The spectrum of eigenmodes traveling 
along a slab can be found from the expression for the 
reflection or transmission coefficients: they have singu- 
larities at the spectral points. For material parameters 
satisfying /i r = — 1 and e r = — 1 the transmission coef- 
ficient equals simply exp(ac?), where a is the decay fac- 
tor of an evanescent mode from the source and d is the 
slab thickness l). Thus, there is only one eigenmode, 
and that solution corresponds to infinitely large k t and 
a. For larger values of a the excitation is closer to the 
resonance with this eigenmode, and the field amplitude 
excited in the slab waveguide is larger. 

In paper p| , however, it has been concluded that if the 
width of the slab is limited, the restoration of fields is 
physically meaningless as it involves infinite energy. To 
clarify the behavior of the evanescent fields in a BW slab 
or a finite width, we study the fields from a source that 
creates only evanescent spectrum in the time domain with 
the finite-difference time-domain (FDTD) method. Time 
domain waveforms at suitably chosen observation points 
and snapshot field distributions are calculated. Our sim- 
ulations show that the amplification indeed occurs, and 
it is due to the surface modes excited on the slab bound- 
aries. 

In any backward-wave medium the negative permit- 
tivity and permeability must be dispersive, and here we 
adopt the Lorentz model to account for the frequency 
dispersion. The expressions for the permittivity and per- 
meability are of the form 



e(u) = eo 1 + 
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Here, uiQ e and ujQ m are the resonant frequencies of the ma- 
terial and r e and T m are the damping factors, which are 
assumed to be zero in numerical simulations. This model 
corresponds to a realization of BW materials as mixtures 
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of conductive spirals or omega particles, as discussed in 
0. In this artificial material both electric and magnetic 
polarizations are due to currents induced on particles of 
only one shape, which provides a possibility to realize the 
same dispersion rule for both material parameters, as in 
©. Note that the medium realized by Smith et. al. is 
built using different principles |7j. Numerical techniques 
appropriate to simulations of a material with the above 
parameters with FDTD can be found in || 0, 0| . The 
discretization scheme used in this paper is described in 
detail in |To| . 

Let the interface between free space and a BW slab lie 
along the x-axis. The excitation plane (a line in our 2D 
cut) is located at a distance d s from the BW slab of thick- 
ness d. In order to be able to study the behavior of the 
evanescent field, we excite a BW slab by a source which 
produces no traveling waves in the direction orthogonal 
to the source plane. We choose the incident electric field 
that depends on the x-coordinate along the interface and 
on the time t in the following manner: 



E(x,t) 



I r(t) cos (k x x — LjQt). 
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The ramp function r(t) increases smoothly from to 1 
over about 50 periods of the cosine function. It is very 
important that the amplitude grows slowly enough; a 
rapid increase in the amplitude does not yield a con- 
stant amplitude on the second slab interface or in the 
image plane, because of a wide frequency spectrum of the 
source, xq is taken to be the x-coordinate in the middle 
of the slab and Xd determines the rate of decay of the 
incident electric field amplitude as measured from xq . loq 
is the center angular frequency of the excitation. Due to 



the finite simulation space, we have given a spatial pro- 
file to the incident field to obtain a decaying amplitude 
near the boundaries of the simulation space. The shape 
of the profile is the normal distribution. The incident 
fields propagate to +£-direction along the interface but 
decay exponentially in the y-direction normal to the slab 
boundaries. From the dispersion relation in free space, 
kx + ky = uj 2 /c 2 we see that by choosing k x > ljq/c 
(k y = 0) in 10} we obtain exponentially decaying fields 
away from the source in the y-direction. In our numerical 
simulations, we use k x = 13.62 m _1 (fco = 9.17 m _1 ). 

The problem space is two-dimensional, with the field 
components H x ,H yi and E z . The peak of the incident 
spectrum is at uiq — 2.75 • 10 9 rad/s, and the parame- 
ters in are the following: ojQ e = u>Q m = 0.55 • 10 9 
rad/s, w*. = u 2 pm = 1.482 • 10 19 (rad/s) 2 , T e = T m = 0. 
With these choices, we obtain e(u) = fi(u) for all uj and 
e(u)/eo = /i(uj)/iio = — 1 at uj = ujq. Absorbing bound- 
ary conditions (ABC) are used to terminate the compu- 
tational domain at the outer boundaries of the lattice. 
For simplicity, we have used Liao's third order ABC 01 > 
although more sophisticated ABC's are available. The 
use of usual ABC's requires a small gap between the 
outer boundary of the computational space and the BW 
material slab, since the Liao's ABC is not applicable to 
dispersive media. The reflection from the slab ends is 
negligible. We excite a slab of thickness 7Ay = 9 cm, 
carrying out the simulations in a 500 x 60 lattice. The 
source plane is located at a distance d s = 3. 5 Ay = 5.25 
cm from the boundary of the slab (tangential magnetic 
field components are defined on the interfaces). Deter- 
mining the precise spatial profile of the source, we take 
xq = 250 Ax and Xd = xq/5 in J3J|. 
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Time step n = 100 



Time step n = 3000 
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FIG. 1: (a) The source fields are visible above the slab, (b) The evanescent fields above the slab create surface waves on 
the boundaries of the slab, (c) The evanescent mode has been amplified by the slab: strong field amplitudes are visible on 
the second interface. The fields created by the source are barely visible above the slab, although they are just reaching the 
maximum amplitude. 
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The problem formulation described above enables us 
to study how the evanescent fields behave in a BW slab 
and how the surface modes are excited. The numerical 
results are presented and discussed next. 

FDTD simulations have been run to see how the elec- 
tric fields behave on the slab boundaries and inside the 
slab. The electric fields created by the source at an early 
stage of the simulation are seen in Figure^a). The snap- 
shot electric field distribution in Figure |T]b) shows that 
the fields have passed through the slab and the field am- 
plitudes on the slab boundaries are high. The electric 
field distribution in Figure^c) is recorded at a later mo- 
ment of time. Figure Q]c) reveals that the electric fields 
are much stronger on the second interface than on the 
first interface. Evidently, the excited evanescent mode 
is amplified by the slab. Notice that the source is just 
reaching its maximum amplitude at n = 3000, although 
the fields created by the source are not visible in Figure 
Uc) due to the scaling. 

Next we record the electric fields as a function of time 
at some observation points. All the observation points 
are on the line x — 250Ax: one in the middle of the slab, 
another on the lower interface and the third in air in the 
image plane, which is located at the distance d — d s from 
the second interface. The electric fields at these points as 
functions of time are shown in Figures and Figure El 
shows that the electric field amplitude is larger on the 
lower boundary of the slab than inside the slab. The in- 
cident electric field amplitude in the source plane is equal 
to unity. The source and the observation point inside the 
slab are located at equal distances from the upper slab 
boundary. Therefore, it is expected that the electric field 
amplitude in the middle of the slab is equal to 1. We 
have observed a smaller amplitude in the middle of the 
slab. The incident electric fields and the fields in the 
image plane are compared in Figure [21 From the the- 
ory, is it expected that the electric field amplitude in the 
image plane equals the incident electric field amplitude. 
The fields in the source plane and in the image plane are 
quite close to each other, see Figure |21 

Let us now compare the results in Figure with ana- 
lytical results. The evanescent mode decays away from 
the source plane as a function of distance y by a factor 
of 
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FIG. 2: The electric field amplitude on the second boundary 
of the slab is larger than that of the source field, and also 
larger than the amplitude in the middle of the slab. This 
result verifies that the evanescent fields are amplified by the 
slab. 



where d s is the distance of the source from the upper 
boundary of the slab. Substituting the parameters we 
obtain T 1.69. In the numerical simulations, the max- 
imum ratio of the field amplitudes on the lower interface 
and on the source plane equals 1.56, being of the same 
order as the analytical result. The fact that the calcu- 
lated value is smaller should have been expected, because 
in the numerical model the incident field amplitude de- 
cays from the slab center, while the estimation is for the 
plane-wave excitation. 



until it hits the upper boundary of the slab. In the slab 
of thickness d, the fields are amplified by a factor of 
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Combining these factors, we obtain the transmission co- 
efficient from the source plane to the second interface as 



(7) 



FIG. 3: The evanescent electric fields in the image plane are 
close to the fields in the source plane. This is also in agreement 
with the theory. 

We have numerically demonstrated that the evanes- 
cent modes are amplified in a lossless frequency dispersive 
backward- wave slab. The field amplitude at the other 
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side of the slab is larger than the source amplitude, and 
the "amplification ratio" agrees well with the analyti- 
cal estimation. The numerically observed electric fields 
in the image plane approximately reconstruct the source 
fields above the slab. Thus, our numerical results sup- 
port the conclusion of paper pj based on the theoretical 
analysis of an infinite lossless backward- wave slab. The 
conclusion of about infinite field energy required for 
the restoration of evanescent fields in finite slabs has no 
theoretical ground (since the slab dimensions are finite 
and the field amplitudes are finite everywhere if all the 
wavenumbers in the spectrum are finite) and is in conflict 
with the results of numerical simulations. 
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